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Abstract 

The relation between thermodynamic phase transitions in classical systems and topological 
changes in their configuration space is discussed for two physical models and contains the first 
exact analytic computation of a topologic invariant (the Euler characteristic) of certain submani- 
folds in the configuration space of two physical models. The models are the mean-field XY model 
and the one-dimensional XY model with nearest-neighbor interactions. The former model under- 
goes a second-order phase transition at a finite critical temperature while the latter has no phase 
transitions. The computation of this topologic invariant is performed within the framework of 
Morse theory. In both models topology changes in configuration space are present as the potential 
energy is varied; however, in the mean-field model there is a particularly "strong" topology change, 
corresponding to a big jump in the Euler characteristic, connected with the phase transition, which 
is absent in the one-dimensional model with no phase transition. The comparison between the two 
models has two major consequences: i) it lends new and strong support to a recently proposed 
topological approach to the study of phase transitions; ii) it allows us to conjecture which particular 
topology changes could entail a phase transition in general. We also discuss a simplified illustrative 
model of the topology changes connected to phase transitions using of two-dimensional surfaces, 
and a possible direct connection between topological invariants and thermodynamic quantities. 
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I. INTRODUCTION 



One can wonder whether the current mathematical description of thermodynamic phase 
transitions (based on the loss of analyticity of thermodynamic observables ^ 0, [J) is the 
ultimate possible one, or whether a reduction to a deeper mathematical level is possible. 

Besides a purely theoretical motivation, there are other reasons for thinking of such a 
possibility. Among the others, we mention the growing experimental evidence that phase 
transitions occur in very small N systems, like nuclear clusters as well as atomic and molec- 
ular clusters, in nano and mesoscopic systems, in polymers and proteins, in very small drops 
of quantum fluids (BEC, superfluids and superconductors). 

Moreover, new mathematical characterizations of thermodynamic phase transitions could 
well be of interest for the treatment of other important topics in statistical physics, as is the 
case of amorphous and disordered systems (like glasses and spin-glasses), or to incorporate 
also first-order phase transitions. A different attempt then is discussed here has been made 



in macroscopic parameter space instead o 



In a number of recent papers 



□ q □ □ n 



in microscopic phase space in Ref. 



10| a proposal has been put forward for 



a new mathematical approach to the study of phase transitions. This applies to physical 
systems described by continuous variables - qi and pi, i = 1, . . . ,N - entering a standard 
Hamiltonian 

1 - 

n = -J2p 2 i+ v (<n>- > (!) 

i=l 

where V(qi, . . . , q^) is the potential energy. The main issue of this new approach is a 
Topological Hypothesis (TH). The content of the TH is that at their deepest level phase 
transitions of a system are due to a change of the topology of suitable submanifolds in its 
configuration space. More precisely, once the microscopic interaction potential V(qi, . . . , g^v) 
is given, the configuration space of the system is automatically foliated into the family 
{E^j^gR of equipotential hypersurfaces independently of any statistical measure we may 
wish to use. Now, from standard statistical mechanical arguments we know that the larger 
the number N of particles, the closer to some are the microstates which significantly 
contribute to the statistical averages of thermodynamic observables. At large N, and at any 
given value of the inverse temperature /3, the effective support of the canonical measure is 
narrowed very closely to a single hypersurface E„ e {g e K.^ | V^(g) = v} C R , with v a 
suitable function of f3. 



Now, the TH consists in assuming that some suitable change of the topology of the {£„}, 
or equivalently, at large N, of the submanifolds M v = {q £ R Ar |^/(g) < v} (the manifolds 
£„ are the boundaries of the M v , i.e., £„ = dM v ), occurring at some v c = v c ((3 c ) (or 
v c — v c {E c )), is the deep origin of the singular behavior of thermodynamic observables at a 
phase transition; (by change of topology we mean that {Y, v } v<Vc are not diffeomorphic 23)] 
to the {H v } v>Vc , or equivalently that {M v } v<Vc are not diffeomorphic to the {M v } v>Vc ). In 
the following of the paper we shall consider the topology changes of the M v , because these 
are naturally investigated using Morse theory. 

The present paper contributes to this new approach to phase transitions with a crucial 
step forward. This consists of an exact analytic treatment of the topology changes in the 
configuration space of the mean-field XY model and of their relation to the thermodynamic 
phase transition, reported in Sec. |HJ In Sec. IIHI we present an analogous treatment of a 
closely related model, the one- dimensional XY model with nearest-neighbors interactions, 
which does not have any thermodynamic phase transition, and we compare the results with 
those obtained for the mean- field case. In Sec. II VI we present and discuss a geometric model 
built in terms of two-dimensional surfaces which should help the intuition in understanding 
some of the aspects of the much more complex topology changes in A-dimensional config- 
uration spaces. Sec. |V| is devoted to some final remarks and to the discussion of possible 
future developments. The paper is completed by Appendix El which is devoted to the proof 
of some crucial estimates used in the main text. 



II. THE MEAN-FIELD XY MODEL 

In this section we give a complete analytical characterization of the topological changes in 
the configuration space of a model with long-range interactions, the mean-field XY model, 
including the one related to the phase transition. This is, to our knowledge, the first com- 
plete analytical characterization of a topological change in configuration space related to 
a t 
in 

energy, we are able - for the first time for a model of physical relevance - to determine 
completely, in a constructive way, the topology and to give analytical estimates of the Betti 



rermodynamic phase transition: a short account of this result has already been given 
I5I ] . For a family of submanifolds of the configuration space defined by the potential 



numbers, i.e., of fundamental topological invariants 



241 ] . This allows us to compute exactly 
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a topological invariant: the Euler characteristic, which can be defined as a combination of 
the Betti numbers jll]. 

The mean- field XY model Q is defined by a Hamiltonian of the class (|TJ) where the 
potential energy function is 

N N 

2..Y 



V ^ = ^mY1 [1 - wa(<Pi - <Pj)] -h^coscpi. (2) 



i,j=l i=l 

Here ipi G [0, 27r] is the rotation angle of the z-th rotator and h is an external field. Defining 
at each site i a classical spin vector Sj = (cos (fi, sin (fj) the model describes a planar (XY) 
Heisenberg system with interactions of equal strength among all the spins. We consider 
only the ferromagnetic case J > 0; for the sake of simplicity, we set J = 1. The equilibrium 
statistical mechanics of this system is exactly described, in the thermodynamic limit, by 
mean-field theory. In the limit h — > 0, the system has a continuous phase transition, with 
classical critical exponents, at T c = 1/2, or e c = 3/4, where e = E/N is the energy per 
particle [16[. 

Defining the magnetization vector m as m = (m x , m y ) , where 

N 

.X 



m x = T7^ c os^ , (3) 



m v N 



1=1 

1 N 

— J^sin^, (4) 

t=i 

the potential energy V can be written as a function of m as follows: 

V((p) = V(m x ,m y ) = y(l - ml- m? y ) - hN m x . (5) 

The range of values of the potential energy per particle, V = V/N , is then 

1 h 2 

-h<V<- + -. (6) 

The configuration space M of the model is an iV-dimensional torus, being parametrized 
by N angles. We want to investigate the topology of the following family of submanifolds 
of M, 

My = V -1 (-oo, v] = {<p G M : V((p) < v} , (7) 

i.e., each M v is the set {ipi\f =l such that the potential energy per particle does not exceed 
a given value v: this is the same as the M v = M E _ K defined above (v has been rescaled 



by jj because we choose V = V/N as a Morse function in order to make the comparison 
of systems with different N easier). As v is increased from — oo to +00, this family covers 
successively the whole manifold M (M v = when v < —h) . 

According to Morse theory |l2j|, topology changes of M v can occur only in correspon- 
dence with critical points of V, i.e., those points where the differential of V vanishes. This 
immediately implies that no topology changes can occur when v > 1/2 + h 2 /2, i.e., all the 
My's with v > 1/2 + h 2 /2 must be diffeomorphic to the whole M, that is, they must be 
A-tori. Moreover, if V is a Morse function (i.e., it has only non-degenerate critical points) 
then topology changes of M v are actually in one-to-one correspondence with critical points 
of V, and they can be characterized completely. At any critical level of V the topology of 
M v changes in a way completely determined by the local properties of the Morse function: a 
/c-handle is attached 2^ 12 1, where k is the index of the critical point, i.e., the number 
of negative eigenvalues of the Hessian matrix of V at this point. Notice that if there are 
m > 1 critical points on the same critical level, with indices k±, . . . , k m , then the topology 
change is made by attaching m disjoint handles H^ kl \ . . . ,H^ km \ This way, by increasing 
v, the topology of the full configuration space M can be constructed sequentially from the 
M v . Knowing the index of all the critical points below a given level v, we can obtain exactly 
the Euler characteristic of the manifolds M v , defined by 

N 

X (M V ) = £)(-l)V*(M w ) , (8) 

where the Morse number fi^ is the number of critical points of V which have index k [12]. 
The Euler characteristic x is a topological invariant (i.e it is not affected by a diffeomorphic 
deformation of M v ): any change in x{M v ) implies a topology change in the M v . It will 
turn out that as long as h > 0, V is indeed a Morse function at least in the interval 
—h < v < l/2 + /i 2 /2, while the maximum value v = l/2 + h 2 /2 may be pathological in that 
it may correspond to a critical level with degenerate critical points. However, as we shall 
show in the following, we shall be able to extend our analysis also to this last critical level. 

Thus, in order to detect and characterize topological changes in M v we have to find the 
critical points and the critical values of V, which means solving the equations 

^=0, i = l,...,N, (9) 

dipt 

and to compute the indices of all the critical points of V, i.e., the number of negative 



eigenvalues of its Hessian 



d 2 V 



N. 



difidifj 

Taking advantage of Eq. (jSJ), we can rewrite the equations Q as 
{m x + h) sin (fi — m y cos (fi — , i — 1, . . . , N 



(10) 



(11) 



As long as (m x + h) and m y are not simultaneously zero (the violation of this condition is 
possible only on the level v = 1/2 + h 2 /2), the solutions of Eqs. (fTTj) are all configurations 
in which the angles are either or tt. In particular, the configuration 



<Pi 







Vz 



(12) 



is the absolute minimum of V, while all the other configurations correspond to a value of 
v which depends only on the number of angles which are equal to ir. If we denote with n n 
this number, we have that the N critical values are: 



1 - ^ (N - 2n^f 



- A (N - 2rO 



(13) 



Inverting this relation yields function of the level value v. 

1 + h 



mt 



-N ± —\\h- 
2 2 



1 1 ' - \ 



(14) 



where int [a] stands for the integer part of a. We can also compute the number C(n n ) of 
critical points having a given n^, which is the number of distinct binary strings of length iV 
having n n occurrences of one of the symbols, which is given by the binomial coefficient 

'N\ N\ 



C(n n ) 



(15) 



We have thus shown that as v changes from its minimum —h (corresponding to n n = 0) 
to | (corresponding to = y) the manifolds M v undergo a sequence of topology changes 
at the N critical values v(n n ) given by Eq. ()13|). We expect that there is another topology 
change located at the last (maximum) critical value, 



1 h 2 

2 + y 



(16) 



However, the above argument does not prove this, since the critical points of V corresponding 
to this critical level may be degenerate. This, because the solutions of the two equations 



in N variables m x = m y = need not to be isolated, so that then on this level, V would 
not be a proper Morse function. Then a critical value v c is still a necessary condition for 
the existence of a topology change, but it is no longer sufficient ^J. However, as already 
argued in Refs. j?i it is just this topology change occurring at v c given in Eq. 
which is related to the thermodynamic phase transition of the mean-field XY model. For, 
the temperature T, the energy per particle e and the average potential energy per particle 
u = (V) obey, in the thermodynamic limit, the equation 

2e = T + 2u(T) ; (17) 

substituting in this equation the values of the critical energy per particle and of the critical 
temperature we get 

u c = u{T c ) = 1/2 . (18) 

as h — > 0, v c — > ~, so that v c = u c . Thus a topology change in the family of manifolds M v 
occurring at this v c , where v c is independent of N, is connected with the phase transition in 
the limit N — > oo, and h — *■ 0, when indeed thermodynamic phase transitions are usually 
defined, at least in the canonical ensemble. 

We still have to prove that a topology change at v c actually exists. To carry out this proof, 
we will use Morse theory to characterize completely all the other topology changes, occurring 
at v < v c : this, together with the knowledge that at v > v c the manifold M v must be an iV- 
torus, will allow us not only to prove that a topology change at v c must actually occur, but 
also to understand in which way it is different from the other topology changes, i.e., those 
occurring at < v < v c . Morse theory allows a complete characterization of the topology 
changes occurring in the M^'s if the indices of the critical points of V are known. In order 
to determine the indices of the critical points (that is the number of negative eigenvalues of 
the Hessian of V at the critical point) we proceed as follows. Since the diagonal elements of 
the Hessian are 

1 1 

Ha = di = — \{m x + h) cos (f { + m y sin ^] - — , (19) 

and the off-diagonal elements are 

H v = ~Jp ( sin fi sin + cos <Pi cos <fj) , (2°) 
one can write the Hessian as the sum of a diagonal matrix D whose nonzero elements are 
Si = — [(m x + h) cosfi + m^sin^j] , i = l,...,N , (21) 



and of a matrix B whose elements are just the Hij given in Eq. (|20p . also for i = j (the 
diagonal elements being —1/N 2 ). Since the ratio between the elements of B and those of 
D is 0(1/N), one would expect at first sight that only the diagonal elements survive when 
N ^> 1, so that the Hessian approaches a diagonal matrix equal to D. However, this is not, 
in principle, necessarily true: one cannot immediately say that at large N the eigenvalues of 
the Hessian are the 8's given in Eq. (|21|) plus a correction vanishing as N — > oo, because the 
number of elements of B is N 2 , so that the contribution of the matrix B to the eigenvalues 
of the Hessian does not, in general, vanish at large N. That nevertheless the argument for 
this crucial point is correct in this special case is shown in Appendix [XJ and is due to the 
particular structure of the matrix B. The latter is of rank one, and one can prove then that 
at the critical points of V the number of negative eigenvalues of H equals the number of 
negative diagonal elements S ±1, so that as ^ > 1 we can conveniently approximate the 
index of the critical points with the number of negative 8's at x, 

index (x) ~ #($ < 0) . (22) 

At a given critical point, with given n r , where the x-component of the magnetization 
vector is 

mx = 1 ~ ~W~ ' ^ 
so that m x > (resp. < 0) if n n < y (resp. > —), the eigenvalues of D are 

Si = m x + h i = 1, . . . , N — n n ; (24a) 

Si = -(m x + h) i = N - n n + 1, . . . ,N. (24b) 

Then, if the external field h is sufficiently small, 

N 

K + h) > if n n < — , 

iV 

{m x + h) < if n w > — , (25) 

so that, denoting with mdex(n n ) the index of a critical point with angles equal to tt, we 
can write 

index(n 7r ) = n n if n v < y , (26a) 

index(n 7r ) = N — n n if > y . (26b) 
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From these equations combined with Eq. (jl5j) one can obtain for the Morse numbers fik, 
i.e., for the number of critical points of index k, as a function of the level v , as long as 
—h < v < 1/2 + h 2 /2 (i.e., excluding the limiting value v = 1/2 + h 2 /2) the following 
expression: 

MW) = (I) ^ ~ " ni ~ ){v)) + e ( N ~ k - n {+) (v))} , (27) 

where B(x) is the Heaviside theta function and n^(v) are the limits of the allowed n^'s for 
a given value of v, i.e., from Eq. ([Ti|). 

iV 



n (±)( w ) = _ 
v 7 2 



I + /; ± i / // 2 — 2 ( u — - 



(28) 



We note that < ni ' < f and f + 1 < nV ] < N, so that Eq. (j2Zj) implies 



N 

fx k (v) = VA;>y, (29) 

i.e., no critical points with index larger than N/2 exist as long as v < v c = 1/2 + h 2 /2. This 
is the crucial observation to prove that a topology change must occur at v c . For, as the Betti 
numbers of a manifold are positive (or zero) numbers, using the Morse inequalities, which 
state that the Morse numbers are upper bounds of the Betti numbers i.e. 

b k <fik for k = 0,...,N, (30) 

we can immediately conclude that as long ast><t> c = l/2 + h 2 /2 

N 

h(M v )=0 Vk> — . (31) 

Thus, as|<f<| + ^- the manifold is only "half" an iV-torus, and since we know that for 
is a (full) iV-torus, whose Betti numbers are 



h(T N )=[f) fc = 0,l,...,JV, (32) 



k 



we conclude that atv = v c = ^ + ^-a topology change must occur, which involves the 
attaching of (^) different /c-handles for each k ranging from f + 1 to JV i.e. a change of 
O(N) of the number of Betti numbers. 

Let us remark that such a topology change not only exists: it is surely a "big" topology 
change, for all of a sudden, "half" an iV-torus becomes a full iV-torus, via the attaching for 
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each different k (ranging from N/2 + 1 to N) of (^) fc-handles. More precisely, a number 
of Betti numbers which is O(N) changes, and changes by amounts which are of the same 
order as their maximum possible values. On the other hand, all the other topology changes 
correspond to the attaching of handles of the same type (index), in fact, as long as v < v c , 
each critical level contains only critical points of the same index, and the index grows with 
v, i.e., if x c and x' c are critical points and V(x' c ) > V(x c ), then index(a^) > index(x c ). The 
potential energy per degree of freedom V is a regular Morse function (or a Morse-Smale 
function jl^) as long as v < v c , but this is no longer true as v > v c ; actually, as we have 
already observed, V could even be no longer a Morse function at all, because the level v c 
might contain degenerate critical points. Nevertheless, as we have shown, this does not 
prevent us from giving a complete analysis of the topology of the M„'s for all the values of 
v, since we can exploit our explicit knowledge of the topology of the M„'s for any v > v c . 

To illustrate what we have described so far, the Morse numbers are shown in Fig. Q as 
a function of k for two values of v, v — |, i.e., an intermediate value between the minimum 
and the maximum of V - shown in Fig. d (a) - and v = | - shown in Fig. E] (b). We see 
that the with < k < y grow regularly as v grows until v c = ^-y-, while all the with 
k > y remain zero, so that also the corresponding Betti numbers must vanish. But at v c 
a dramatic event occurs, because for all the values of v > v c the Betti numbres must be 
those of an iV-torus, which are reported for comparison in Fig. d (c). A sudden transition 
from the situation depicted in Fig. [T] (b) to that of Fig. [T] (c) occurs at v c , i.e., y Betti 
numbers simultaneoulsy become nonzero. 

Now we are in a position to explain how these topological transitions can be described 
by topological invariants. The optimal situation would be the possibility of computing all 
the N + 1 Betti numbers of the manifolds M v as a function of v: unfortunately, we are only 
able to set an upper limit to them, using the Morse inequalities flUOj) . Nevertheless, we can 
use Eqs. (jHJ, (jHj) and (|27j) to compute the Euler characteristic of the manifolds M v , since 
we only need the for that. It turns out then that x jumps from positive to negative 
values, so that it is easier to look at In Figure[21 log(\x\(M v )) / N is plotted as a function 
of v for various values of N ranging from 50 to 800. The "big" topology change occurring 
at the maximum value of V, which corresponds in the thermodynamic limit to the phase 
transition, implies a discontinuous jump of \x\ from a big value to zero. 
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III. THE ONE-DIMENSIONAL XY MODEL 



We consider now one example where there are topological changes very similar to the 
ones of the mean- field XY model but no phase transition occurs, i.e., the one-dimensional 
XY model with nearest-neighbor interactions, whose Hamiltonian is of the class (JTJ with 
interaction potential 



In this case the configuration space M is still an A-torus, and using again the interaction 
energy per degree of freedom V = V/N as a Morse function, we can prove that also here 
there are many topological changes in the submanifolds M v as v is varied from its minimum 
to its maximum value. 

The critical points are again those where the (pi s are equal either to or to ir. However, 
at variance with the mean-field XY model, the critical values are now determined by the 
number of domain walls, rid, i-e., the number of boundaries between connected regions on 
the chain where the angles are all equal ("islands" of 7r's and "islands" of O's). The number 
of 7r's leads only to a correction 0(h) to the critical value of v, which is given by 



Since rid £ [0, A — 1] (with free boundary conditions, rid = 0,1, N — 1, while with periodic 
boundary conditions rid is still bounded by and N — l, but can only be even), the critical 
values lie in the same interval as in the case of the mean-field XY model. But now the 
maximum critical value, instead of corresponding to a huge number of critical points, which 
rapidly grows with N, corresponds to only two configurations with N — l domain walls, 
which are y?2fc = 0, y?2fc+i = 7T, with k = 1, . . . , N/2, and the reversed one. 

The number of critical points with rid domain walls is therefore (assuming free boundary 
conditions) 



We can compute the index of the critical points also in this case (see Appendix IA2I for 
details). It turns out that 



N N 




(33) 




(34) 




(35) 




(36) 
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so that 



ix k {n d ) = 2 6(n d - k) 




(37) 



It is evident then that any topology change here corresponds to the attaching of handles 
of the same type. However, no "big" change like the one at v c in the case of the mean- 
field model exists, although V is a Morse-Smale function on the whole manifold M. To 
illustrate this, we plot in Fig. El the values of the Morse indices fik as a function of k, as we 
have already done for the mean-field XY model in Fig. ^ Comparing Fig. E] with Fig. 
we see that while in the mean-field model there is a critical value where y Betti numbers 
become simultaneously nonzero, i.e., there exists a topology change which corresponds to 
the simultaneous attaching of handles of y different types, while here in the one-dimensional 
nearest-neighbor model nothing like that happens. 

Also in this case, using Eqs. (|3*5|) and we can compute the Euler characteristic of 
the submanifolds M„: 



The Euler characteristic for the one- dimensional nearest-neighbor case is shown in Fig. 0] 
Comparing this figure with Fig. El we see that there is here no jump in the Euler character- 
istic. 

Before discussing the relevance that these results may have for the general problem of 
the relation between topology and phase transitions, it is illustrative to consider two ab- 
stract simplified models of topological transitions that occurred in the two models we have 
considered so far. 

IV. TWO-DIMENSIONAL MODEL OF TOPOLOGY CHANGES 

A two-dimensional model of the topological transition occurring in the configuration space 
of the physical models we have discussed, which could perhaps help the intuition, can be 
built as follows. Let us consider a two-dimensional torus T, and place it on a plane: let /i max 
be the maximum height of the surface above the plane (see Fig. El (a)). Then deform the 



(38) 



where, due to Eq. lJ33j) . 



n d (v) = 2Nv + 0(h) . 



(39) 
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torus (by means of a diffeomorphism) until the upper end of the hole is at height h T 
obtaining the surface M shown in Fig. |3] (b) . It is apparent that e can be made as small as 
we want. 

Let us now consider the height function 7i above the plane as a Morse function, so that 
the manifolds 

T h = H~ l (-oo, h] = {x G T : H{x) < h} (40) 
M h = H~\-oo, h} = {xeM : H{x) < h} (41) 

are defined. As h varies from its minimum to its maximum values (h = and h — 3, 
respectively, in Fig. |SJ), the manifolds 1\ and Mh cover the whole torus; as long as h is lower 
than the top of the hole (h = 2 in Fig. EJ), both and Mh are "half-tori", but then the 
Th become gradually a full torus, while the Mh jump abruptly from a half-torus to a full 
torus as h is changed by e. Identifying the height function with the potential energy, the 
case of the clearly reminds of the behavior of the one- dimensional XY model with no 
phase transition, while the case of the Mh seems close to what happens in the mean-field 
XY model, and the "jump" from the half-torus to the full torus is similar to the topology 
change which is connected to the phase transition. 

The analogy with the mean-field models becomes even clearer if, instead of a torus, we 



consider a compact surface of genus |26j g S> 1, i.e., with many holes, and deform the surface 
with a diffeomorphism until the upper end of all the holes is at height h m3jX — e, as shown 
in Fig. H3 Again, e can be made as small as we want. Let us denote by M^ the deformed 
surface. 



22 



The Betti numbers of the surface M^ 9 ' are 

b = b 2 = 1 , h = 2g, (42) 

and the Euler characteristic is 

x(MW)= X (MW) = 2-2 Sl (43) 

i.e., a big negative number. Let us now consider, as in the previous case of the torus, the 
height function TC above the plane as a Morse function. The manifolds 

M {g) = W -i(_ 0Gj h] = {xe M (9) : H{x) < h} (44) 
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will be topologically very different from the whole M^ 9 > as long as h < /i max -e but sufficiently 
large that all the critical levels corresponding to the bottoms of all the holes have already 
been crossed: in fact, their Betti numbers will be 

b (Mi 9) ) = 1, h(Mi 9) ) = g, b 2 (Mi 9) ) = , (45) 

and the Euler characteristic will be 

X (M { h 9) ) = l-g. (46) 

Then, by changing the value of the height by an amount e as small as one wants, one changes 
the Betti number b\ from g to 2g, the b<i from to 1 and x from 1 — g to 2 — 2g. This 
is a topological change which involves a change of 0(d) Betti numbers (d is the dimension 
of the manifold); moreover, the size of the change is of the order of the value of the Betti 
numbers. This topology change involves also a change of the Euler characteristic x which is 
again of the same order as its value. Identifying again the height function with the potential 
energy, we see that this is just what happens in the case of the M v of the mean- field XY 
model, although there the dimension of the manifolds is N and very large, while in this 
low-dimensional analogy it is only d = 2. The behavior of )| as a function of h is 

plotted in Fig. We see that the behavior of |%| is indeed very similar to the case of the 
mean-field XY model, the only big difference being that in the latter \x\ jumps to zero while 
here it jumps to a nonzero value. However this difference is due to the fact that the Euler 
characteristic of a torus is zero while that of a surface of genus g is 2 — 2g. 



V. CONCLUDING REMARKS AND OPEN QUESTIONS 

We conclude with some comments and remarks, and with a discussion of the open prob- 
lems and of the future perspectives opened by the results we have presented. 



A. Topology and thermodynamic functions 

The results we have reported in the present paper concerning the mean-field and the 
one-dimensional XY models, together with a theorem [l^J which states the necessity of 
topology changes in the M„'s for the existence of phase transitions (at least for a certain 
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class of short-range interactions), lend strong support to the TH. However, one could wonder 
whether a direct relation between topology and thermodynamic quantities, like entropy or 
temperature, exists. 

A result - although approximate - which shows that there could be an explicit contribu- 
tion of topological quantities to the temperature was indeed found in a previous paper [10]: 
it relates the inverse temperature (3 = ^ to the Betti numbers 6, of the constant-energy h 



ny- 





persurface in phase space, Eg. We stress the fact that such a relation - Eq. (43) in Ref. 
- does not allow to express (3 in terms of topological invariants alone, because also other 
contributions, of a non-topological nature, exist: nonetheless a topological contribution ex- 
pressed as a function of the Betti numbers of Eg is explicitly present. A similar relation can 
be obtained also for the en trop y: the derivation only involves a slight modification of the 

lpl. Without entering into technical details, we simply observe 



reasoning of Sec. V of Ref. 



that, under the same approximations as made in Ref. Q], the "topological contribution" 
t(E) to the entropy per degree of freedom s = -j| can be written, up to constant terms, as 

1 2JV-1 

r(jE;) ^2iv^T log S 6i(S£) ' (47) 

i=0 

an expression identical to the corresponding expression for f3 in Q], apart from the presence 
of a logarithm here. We remark that, as in the case of f3, there are also other contributions, of 
a non-topological nature, to the entropy. However, it is interesting to look at this topological 
r-contribution alone to see whether it carries the relevant information about the way the 
topological transition of the energy hypersurface in phase space triggers the phase transition, 
i.e., the nonanalyticity of the thermodynamic quantities. Therefore we can try to compute 
t for the models studied in the present paper. 

First of all, we have to express r(E) in terms of the topological invariants of the manifolds 
M v instead of in terms of the Betti numbers of the energy hypersurfaces E# [Eq. (|47|) ]. To 
do that, we resort to the fact that - at large N and for systems with standard kinetic energy 
K - the volume measure of E B concentrates onj27| §>^~j 1/2 x M^ V ), so that E E can be 
approximated by this product manifold. Using then the Kunneth formula Q] for the Betti 
numbers of a product manifold A x B, i.e., 

b t (AxB)= b M)h(B) (48) 

j+k=i 
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which, when applied to 1 



x M(v>, gives bi(E E ) = 26 i (M«) for i = 1, . . . , JV - 1, and 
bj(T,E) = bj(M v ) for j = 0, N, since all the Betti numbers of a hypersphere vanish except b 
and b N which are equal to 1. This allows us to write 

N-l 



t(v) 



log 

2N & 



b (M v ) + 2j2HM v )+b N (M v ) 



i=l 



(49) 



where we have also replaced 2N — 1 with 2N because we assume that N is large. Since we 
cannot compute the Betti numbers exactly, we cannot evaluate this expression; nonetheless, 
we can estimate it by using the Morse numbers as approximations of the Betti numbers, i.e., 
by putting 

h{M v ) ~ m{M v ) , 2 = 0,..., N. (50) 

In general, with the exception of perfect Morse functions (i.e. those for which bi = //j), the 
Morse numbers are only upper bounds of the Betti numbers. However, looking at Figs. ^ 
(b)-(c) (for the mean-field XY model) and El (b)-(c) (for the one- dimensional XY model), 
we can observe that the nonvanishing Morse numbers of the manifolds M v for values of v 
very close to t> max are very close to the Betti numbers of the M„'s for v = v max (which 
are iV-tori), and this suggests that for the models studied in the present paper the Morse 
numbers \i{ could well be good approximations to the Betti numbers 6j. 
We can then compute the quantity 

N-l 

fio(M v ) + 2 fH(M v ) + Hn(M v 



T(V 



2N 



log 



i=l 



(51) 



for the two models studied in the paper. We note that 



N 



J2^(M V )=N C (M V ) , 



(52) 



i=0 



where N C (M V ) is the total number of critical points of the function V in the manifold M v . 
Therefore, when N gets large, and if no and fi^ are not much larger than the other Betti 
numbers (which is true in our case, see Fig. QJ, we can write approximately 

1 



T(V 



N 



\ogN c (v) 



(53) 



apart from an additive constant. The behavior of this quantity is plotted as a function of 
v in Fig. |H]for the mean-field XY model and in Fig. |§]for the one-dimensional XY model. 
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First, we note that in both cases the "topological contribution" f to the entropy behaves 
qualitatively as expected for the configurational entropy, i.e., it grows monotonically up to 
the maximum value of v, after which it remains constant. Moreover, we see that, in the case 
of the mean- field model (Fig. |HJ), the topology change at v c related to the phase transition 
corresponds to a discontinuity in the slope of f(v), which thus seems to be both the precursor 
and the source of the nonanalyticity of the entropy at the phase transition point as iV — > oo. 
In the case of the one-dimensional XY model, where no phase transition is present, unlike 
the mean-field case, the curve is smooth for any v, consistent with the fact that also the 
entropy of the system is smooth. 

Therefore the contribution of the Betti numbers alone to the thermodynamic quantity s 
appears to yield the phase transition point correctly; an indication again of the topological 
mechanism of the phase transition. 

We emphasize that the discontinuity in Fig. |21 at v = v c corresponds to that of Fig. |H] at 
the same value of v. Even more, the many small jumps in the Euler characteristic occurring 
in Fig. El and corresponding to the topological changes which occur at v < v c , are smoothed 
out in Fig. |H1 where the topological contribution to the entropy, f, is reported. This is due 
to the fact that while the Euler characteristic is the alternating sum of the Morse numbers 
fii, the r is the sum of them: this is a further indication that to yield a phase transition, 
i.e., a discontinuity in a derivative of s, a "strong" topology change would be needed, such 
as to affect the variation of the sum of all the Betti numbers as a function of v. 

However, to what extent these results are of general validity remains an interesting open 
question. 

B. When does a topology change entail a phase transition? 

Another fundamental point which still remains open in the topological approach to phase 
transitions is the question of which are the sufficient conditions for the topology changes 
of the manifolds M v to entail a phase transition. Topology changes appear to be rather 
common, and most of them are not connected to phase transitions. 

In our previous papers, clear evidence - albeit only numeric - was found that phase transi- 
tions would correspond to very significant transitions in the way the topology changes |8j as a 
function of v. Moreover, in Sec. IV Al we have discussed a possible mechanism through which 
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a topological change could trigger a change in the properties of thermodynamic functions, 
resulting in a phase transition. 

On the basis of our results obtained before and here, we put forward the conjecture that 
what we have observed in the case of XY models may well have general validity: a topology 
change in the submanifolds M v might entail a phase transition if it involves the simultaneous 
attachment of handles of O(N) different types on the same critical level. However, we must 
be aware that this could well be specific to the class of models studied here: to what extent 
this conjecture might have general validity remains an open question. 



C. The role of the external field h 



Studying the topology changes in the configuration space of XY models, mean-field 
as well as one-dimensional, we have considered the presence of an external field h 7^ 
which explicitly breaks the 0(2)-invariance of the potential energy, and then, discussing the 
connection with phase transitions, we have considered the case in which h tends to zero. We 
did that for the sake of simplicity, for, if we set h = from the outset, the potential energy 
per degree of freedom V is not, rigorously speaking, a Morse function, because its 0(2)- 
invariance entails the presence of a zero eigenvalue in its Hessian. When h = 0, the critical 
points of V are not isolated, but form one-dimensional manifolds (topologically equivalent 
to circles) which are left unchanged by the action of the 0(2) continuous symmetry group 
so that the critical points become in this case critical manifolds. However, in the case of 
the mean-field XY model, as far as the presence and the nature of topology changes are 
concerned, studying the case with h = from the outset we find exactly the same behavior 
as in the case we have discussed in this paper, i.e., as long as v < v c only handles of the 
same type are attached, while at v — v c handles of N/2 different types are attached, the 
only difference between the two cases being that when h = the handles are not attached 



at isolated points, but rather to the entire critical manifold [12j,[13(. However, putting h = 
from the beginning makes the computation of the Euler characteristic x(M v ) via the Morse 
numbers much more difficult, because now one has to take into account the contributions 
to x comulg ho m the Betti n^bers of the critica, — , S ee Re, for the deta*). 

In the case of the one-dimensional nearest-neighbor XY model, where no phase transition 
is present, the use of h = from the outset implies a further complication, i.e., that the 
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critical points consist not only of the configurations made of O's and 7r's, but also of spin 
waves, that is configurations 

<Pj = ip e ikj , (54) 

with wavenumbers k depending on boundary conditions. 

For all these reasons we preferred to force the potential energy to be a Morse function via 
the explicit breaking of the 0(2) symmetry using an external field h ^ 0. Incidentally, we 
notice that in the mean-field XY model, as long as h ^ 0, the topology changes which do not 
correspond to any phase transition (i.e, those occurring at v < v c ) occur at a number of values 
of v which grows with N, and these values become closer and closer as N grows, eventually 
filling the whole interval [0, |] as N — > oo. To the contrary, the value v c , which corresponds 
to the "big" topology change connected to the phase transition, remains separated from the 
others by an amount 0(h 2 ) also in the thermodynamic limit, and tends to \ only when h — > 0. 
This is reminiscent of a similar fact occurring in statistical mechanics, where one observes a 
spontaneous symmetry breaking, signalled e.g. by the onset of a finite magnetization even 
at zero external field, if one assumes the presence of an external field and then lets it tend 
to zero only after the thermodynamic limit is taken. 

D. Transitional phenomena in finite systems and other future developments 

(a) If it would be possible to establish a one-to-one correspondence between a particular 
class of topology changes in configuration space and the usual thermodynamic phase tran- 
sitions defined in the thermodynamic limit, then, as a by-product, one would have available 
also a natural definition of phase transitions for finite systems. In fact, the topology changes 
are defined at any N, so that one would call "phase transitions at finite N" those topology 
changes in configuration space which become phase transitions in the usual, statistical- 
mechanical sense at infinite N. This would be an interesting consequence of the present 
topological approach to phase transitions, since it would circumvent the basic problem of 
any clear definition of phase transitions at finite N, by showing that all phase transitions 
have their basic origin in a topological change which may occur also at finite N, but entail 
true mathematical singularities in the thermodynamic functions only at infinite N. 

(b) Finite systems do not seem to be the only active field of statistical physics where the 
topological approach might prove useful. Another promising field of application is that of 
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glasses or, in general, of disordered systems. It is now quite clear that many of the puzzling 
properties of glasses are encoded in their "energy landscape" [l^i i.e., in the structure of 
valleys and saddles of the potential energy function: but this is directly connected to the 
structure of the submanifolds M v and of configuration space, and in fact topological 
concepts start to emerge in some recent papers on glasses [3]. 

(c) Finally, we notice that at present only systems undergoing second-order phase tran- 
sitions have been studied via the topological approach: a natural question which arises is 
then if also first-order transitions can be explained topologically. Work is in progress in this 
direction, and some preliminary results indicate that also discontinuous phase transitions 
can be connected to topology changes of the submanifolds M v , and that a signature of them 
can be found in their Euler characteristic. 
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APPENDIX A: ESTIMATE OF THE INDEX OF THE CRITICAL POINTS 

In this Appendix we want to discuss some details related to the estimates of the indices 
of the critical points we have used in Sec. Eland Sec. 11111 



1. Mean-field XY model 

Here we want to prove the crucial estimate (|22|) . which we used in Sec. |H]to compute the 
index of the critical points for the mean-field XY model. 

We recall that we want to compute the the number of negative eigenvalues of the Hessian 
matrix of the function V, i.e., of the matrix H whose elements are 

d 2 V 



Hij = ^r^T i,j = l,...,N, (Al) 
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where V = V/N and V is the potential energy of the mean- field XY model defined in 
Eq. (J2J). The diagonal elements of this matrix are 

1 1 

Ha = di = — [(m x + h) cos <pi + m y sin - — , ( A2) 

and the off-diagonal ones are 

Hij = -J-^ ( sin Pi sin ¥j + cos fi cos fj) ■ ( A3 ) 

At the critical points of V, the angles are either or tt, so that the sines are all zero and the 
cosines are ±1. Moreover, since we are interested only in the signs of the eigenvalues of H 
and not in their absolute values, we multiply H by N in order to get rid of the 1/N factor 
in front of it. We can then write the matrix H (multiplied by N) as 



where D is a diagonal matrix, 



whose elements 5 are 



H = D + B (A4) 



D = diag^) , (A5) 



Si = (m x + h) cos ifi , (A6) 



where the ipi's {i — 1, . . . , N) are computed at the critical point, and the elements of B can 
be written in terms of a vector o whose N elements are either 1 or —1: 



hj = --Jj a i a i ' ( A7 ) 



where 

at = +1(-1) if <pi = 0(tt) . (A8) 
This, since when the angles are either or it, the sines in Eq. ()20j) vanish, so that then 

NH i:j = -— cos ifi cos ifj , (A9) 

and 

COSf^j = CTj . (A10) 

Having fixed the notation, our goal is to show that, at least when iV is large, the number 
of negative eigenvalues of the full matrix H, i.e., the index of the critical point, can be 
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conveniently approximated by the number of negative eigenvalues of D, that is, by the 
number of negative 8's. To do that, we proceed in two steps: (i) we show that the matrix 
B is of rank one (which implies that B has N — 1 zero eigenvalues and only one nonzero 
eigenvalue), and (ii) we adapt a theorem due to Wilkinson 17] to this case, thus proving 
our assertion. 

As to step (i), let us consider for example a case with N = 3, and the critical point 
corresponding to, say, (<pi, <^ 2 , (fs) = (vr, 0, 0). The vector a is then 

K <7 2 ,<7 3 ) = (-1,1,1) • (AH) 

Using Eq. (jA7j) . the matrix B is 

/ 1 -1 -A 

. (A12) 

(-11 1 J 

We see that the second row is equal to the first multiplied by —1, and the same holds for the 
third row. This is true for any N, and is a consequence of Eq. (J A 7)1 : any row of the matrix 
B is equal to another row multiplied by either +1 or —1. This means that N — 1 rows are 
not linearly independent and that the rank of the matrix is one. 

We have then proved that our Hessian matrix H is the sum of a diagonal matrix and of 
a matrix of rank one. 

Let us now pass to step (ii). First, we recall a theorem of Wilkinson found in Ref. jl7|: 
Theorem A.l (Wilkinson) Let A and B be N x N real symmetric matrices and let 

C=A+B . 

Let 7j, aj and Pi (i = 1,2, . . . , N) be the (real) eigenvalues of C , A and B, respectively, 
arranged in non-increasing order, i.e., 

7i > 72 > • • • > In ; 

oc\ > «2 > ... > «jv ; 

A > fa > ■■■ > Pn- 

Then 

lr+s-i < a r + p s Vr + s-l<A. (A13) 
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Notice that we can also write 

A = C+(-B), (A14) 
and since the eigenvalues of —B arranged in non-increasing order are 

we can write 

«r+ S -l < 7r - Pn-s+1 ■ (A15) 

Now, we are interested in a special case, i.e., the case in which the matrix B is of rank 
one (has only one nonzero eigenvalue). What does Wilkinson's theorem say when applied 
to such a special case? We consider the two possible cases, namely: 

(a) the nonzero eigenvalue is negative: 

Pi = 0foTi = l,2,...,N-l,l3 N = -Q; 

(b) the nonzero eigenvalue is positive: 

P N = g, Pi = for i = 2,3,. . ., N. 

Case (a). Choosing s = 1 in Eq. (jA13|) we get ft s = 0, so that 

lr <a r r = l,...,N, (A16) 
while choosing s = 2 in Eq. (|A15|) we get — PN — s + 1 = again, whence 

«r+i < 7r r = 1, . . . , N - 1 . (A17) 
Combining Eqs. (jA16|) and (jA17|) we obtain 

Or+i < 7r < Or r = 1, . . . , N - 1 ; (A18) 

lN<a N . (A19) 

We have thus shown that all the eigenvalues of C (except for the smallest one) are bounded 
between two successive eigenvalues of A. As to 77V, we can only say that it is smaller than 
(or equal to) the smallest eigenvalue of A. 
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Case (b). Choosing s = 2 in Eq. (|A13|) we get (3 S = 0, so that 

Tr+i < r = 1, . . . , N - 1 , (A20) 

while choosing s = 1 in Eq. (jA15j) we obtain 

«r<7r r = l,...,iV. (A21) 
Combining Eqs. (|A20|) and (jA21j) we obtain 

a r < lr < a r _i r — 2, . . . , N ; (A22) 

7i > ai . (A23) 

We have thus shown again that all the eigenvalues of C (except, in this case, for the largest 
one) are bounded between two successive eigenvalues of A. As to 71, we can only say that 
it is larger than (or equal to) the smallest eigenvalue of A, a±. 

Let us now apply these results to our problem, i.e., to the computation of the number 
of negative eigenvalues of the matrix H in Eq. (jA4|) . Denoting by rji its eigenvalues, by 5i 
those of D and by those of B, we have 

Pi = Vz^l,iV, (A24) 

and either 

fix = 0,P N = -g, 

or 

Pi = g, (3 N = . 

At a given critical point, with n n angles equal to tt, the eigenvalues of D are (for the moment 
we do not order them) 

5i = m x + h i = 1, . . . , N — n n ; (A25) 

S i = -(m x + h) i = N - n n + 1, . . . , N. (A26) 
The x-component of the magnetization vector is given by 

m x = 1 - ^ (A27) 
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so that 

N 

m x > if n„ < — , (A28) 
N 

m x < if TV > — . (A29) 

Then, if the external field h is sufficiently small: 
if < N/2, then 

5i = m x + h > i = 1, . . . , N - n w , (A30) 



Si = -(m x + h) < i = N - n n + 1, . . . , N, (A31) 

i.e., there are N — n n positive and negative 5's; 
else if < iV/2, then 

5- = -(m x + h) > i = 1, . . . , n„. , (A32) 



^ = m x + /i< i = TV + 1, . . . , N , (A33) 

i.e., there are positive and N — n n negative 8's. 

Now we claim that, at least as gets large, we can estimate the number of negative 77's, 
i.e., the index of the critrical point, by saying that it is equal to the number of negative 5's. 
More precisely, we claim that the error of our estimate is not larger than 1, i.e., 

mdex(H) = #(77 < 0) = #(<J < 0) ± 1 . (A34) 

and as gets large this error becomes obviously negligible. To prove this statement, let 
us consider the case in which n n < N/2. We observe that we do not know whether we are 
in case (a) or in case (b), i.e., we do not know if the matrix B has a negative or a positive 
eigenvalue. But we can, anyway, try one of the two cases, say (a). Using Eqs. (|A18J) and 
(IA19|) we can then say that 

S r +i <T] r <5 r <0 r = N-n 7T + l,...,N-l, (A35) 

(note that these are n n — 1 equations), and that 

r] N <5 N <0. (A36) 
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Thus we conclude that the number of negative 77's is just equal to that of negative S's, i.e., 
n w . If we guessed correctly the sign of the nonzero eigenvalue of B, then, our estimate is 
exact. But in case we guessed it wrong, i.e., if we were in case (6) and not (a), by using Eqs. 
flA18|) and (jA19|) . we would have overestimated the number of negative 77's by 1. Conversely, 
if we had used the equations of case (b) in a situation which belonged to case (a) we would 
have underestimated the index by 1. So, we conclude that the error of our estimate is always 
±1. 



2. One-dimensional XY model 

Here we want to discuss the details of the result reported in Eq. (JHEj) . i.e., that in the 
case of the one- dimensional XY model with nearest-neighbor interactions the index of the 
critical points equals the number of "domain walls" in the configuration. 

First of all, let us notice that in the present case the Hessian matrix H is tridiagonal, i.e., 
it can be written as 

/ a x (3 X ... \ 
... 



H 



Pi a 2 P2 

••• '•• 







. . . (3 N - 2 a N -i /3jv-i 
(3 N -! a N J 



(A37) 



\ ... 

where, assuming free boundary conditions, 

a.\ = cos(<£> 2 — (fx) + h cos(v?i) ; (A38a) 
cxj = cos(y? m - ifi) + cos(y9i - <^_i) + h cos(y?i) , % = 2, . . . , N - 1 ; (A38b) 

OL TV = COs((y9/v — <^AT-l) + /iCOs(v9/v) , (A38c) 

and 

fa = - cos(<p i+1 - <fi) , i = l,...,N-l. (A39) 
Since at critical points tpi = or it, we have that for any i and for any critical point 



(A40) 
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while the diagonal elements «j are 



at = l±h; (A41a) 
Oi = 2±h, i = 2,...,N- 1 ; (A41b) 
a N = 1 ± h , (A41c) 

if there are no domain walls, i.e., if = 0, while they can assume also the values ±/t and 
—2 ± h (—1 ±/iifz = lorz = A) if 7^ 0, i.e., if there are domain walls. 

Let us now prove that 7^ is a necessary condition for the presence of negative 
eigenvalues of the Hessian, i.e., for a nonvanishing index of a critical point. To do that, we 



recall a theorem due to Gershgorin (see e.g. Ref. 
symmetric matrix, can be stated as follows: 



i7j; 



which, in the simple case of a real 



Theorem A. 2 (Gershgorin) Let A be a real n x n symmetric matrix whose elements are 
a^, and let 



7~j ^ ^ I CLij I , i 1 j • • • , Tl 



the eigenvalues of A lie in the intervals 

Xi = {x G K : \x- a u \ < r^} ; 

if m of the Xi form a disjoint set, then precisely m eigenvalues (counted with their multi- 
plicity) lie in it. 

In our case, due to Eq. (jA40|) . at any critical point we have 

ri = r N = 1; (A42a) 
n = 2 , i = 2, . . . , N - 1 , (A42b) 

so that, if Ha = and h — > 0, Eqs. (|A41|) and Gershgorin's theorem imply that all the 
eigenvalues lie in the interval \x — 2| < 2, so that there are no negative eigenvalues and the 
index is zero. On the other hand, if 7^ and h — *■ 0, then the intervals Xi are either 
\x\ < 2, or \x + 2| < 2, hence the eigenvalues lie in the interval (—4,2), so that the index 
can be nonvanishing. However, Gershgorin's theorem is useless to compute the number 
of negative eigenvalues, because the intervals Xi overlap each other, thus the eigenvalues 
cannot be localized more strictly. 
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Anyway, the fact that the Hessian is tridiagonal allows us to compute directly its char- 
acteristic polynomial det(H — XI), whose roots Ai, . . . , Xn are the eigenvalues, by means of 
a recurrence formula. Let 



Po(A) = 1 ; 
Pi(A) = «i - A ; 

p fc (A) = (a k - X)p k -i(X) - /3fc_ 1 p fc - 2 (A) 



(A43a) 
(A43b) 
(A43c) 



then, since 



p k (X) = det 



( a x -X fix 
Pi a 2 -X p 2 

o ••. ••. ■•. 



\ 







k = 2, . . . , N , (A44) 



... p k _ 2 afc-1 - A p k -i 
\ ... p k - x a k -Xj 

the characteristic polynomial of H is given by Pn{X). Since at the critical points all the /3's 
are ±1 - see Eqs. (|A40|) . we have that 



Po(A) = 1 ; 
Pi(A) = ai - A ; 

p fc (A) = - A)pfc-i(A) -Pk-2(X) 



(A45a) 
(A45b) 
(A45c) 



so that the characteristic polynomial Pjv(A) depends only on the a's. Moreover, the following 
theorem holds (see e.g. Ref. |2l|): 



Theorem A. 3 Let H be a tridiagonal symmetric matrix defined as in Eq. \A3Ty . Define 
the sequence 

{p (X), Pl (X),...,p N (X)} (A46) 

as in Eq. \AJj.Sfy ; then the number of sign changes in the sequence (with the rule that if 
Pi(X) = then it has the opposite sign of pi^i(X)) equals the number of eigenvalues of H 
which are less than or equal to X. 



Then the number n c of sign changes in the sequence 

{p (0),pi(0),...,^(0)} 



(A47) 
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equals the number of negative eigenvalues, i.e., the index of the critical point because no 
eigenvalues are zero. If one puts h = 0, then there is one eigenvalues which becomes zero 
at any critical point, so that the index equals n c — 1, but in this case one easily sees by 
direct computation (which can be performed exactly on a computer at any N because in 
this case the a's are integer numbers) that rid = n c — 1, so that one finds the result reported 
in Eq. flUJ). 
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FIG. 1: Mean-field XY model, (a) Histogram oflog(nk(M v ))/N as a function of k for v = 1/4; (6) 
Histogram of log(fik(M v )) / N as a function of k for v = 1/2. In both cases N = 50 and /i = 0.01. (c) 
For comparison, histogram of log(6fc(T iV ))/A r as a function of k for an iV-torus T N , with N = 50, 
which is the lower bound of log(^(M„)) / 'N for any v > v c . 
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FIG. 2: Mean-field XY model. Plot of log(|x|(M„))/iV as a function of v. N = 50,200,800 (from 
bottom to top) and h = 0.01; v c = 0.5 + 0{h 2 ). 
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FIG. 3: The same as Fig. {I\ioT the one-dimensional XY model with nearest-neighbor interactions, 
(a) Histogram of \og(nk(M v )) / N as a function of k for v = 1/4; (6) Histogram of log(/j 1 f i .(M v )) / N 
as a function of k for v = 1/2. In both cases N = 50 and h = 0.01. (c) For comparison, histogram 
of log(bk(T N ))/N as a function of k for an iV-torus T^, with N = 50. 
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FIG. 4: Plot of log(|x| (M„))/iV for the one-dimensional XY model with nearest-neighbor interac- 
tions as a function of v. N = 50, 200, 800 (from bottom to top). 




FIG. 5: (a) A torus T and its height function. Here /i max = 3. (b) A deformation M of such a 
torus as explained in the text. 
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FIG. 6: (a) A surface of high genus (here g = 9). (6) A deformation of such a surface as 

explained in the text. 




FIG. 7: Plot of the logarithm of the absolute value of the Euler characteristic of the submanifolds 
of given height of a surface M^ 9 ' like the one depicted in Fig. |B] (b), with g = 50 and h m3X = 3, 
as a function of the height h. The small jumps are the topology transitions corresponding to the 
crossing of the bottoms of the holes: the last big jump is the one occurring at h max , when e — > 0. 
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FIG. 8: Mean-field XY model. Plot of log(N c )/N as a function of v. N = 50,200,800 (from top 
to bottom) and h = 0.01; v c = 0.5 + 0(h 2 ). 




FIG. 9: One-dimensional XY model with nearest-neighbor interactions. Plot of log(N c )/N as a 
function of v. N = 50,200,800 (from top to bottom) and h = 0.01. 
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